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1.  Introduction 

Models  frequently  used  for  the  hydraulic  machines  are  composed  of  a  disc  and  two  perpen¬ 
dicular  radial  plates  in  the  flow  produced  either  by  a  linear  source  (Fig.  la)  or  by  a  linear  pit.  Or, 
neglecting,  in  the  first  approximation,  the  interaction  between  the  disc  and  two  plates,  in  this  case 
two  boundary  layers  can  be  considered  separately  [8]:  the  planar  boundary  layer  between  two  di¬ 
vergent  (or  convergent)  plates  and  the  axisymmetric  boundary  layer  on  the  disc  (in  the  presence  of 
a  funnel  above  the  disc  (Fig.  lb)  to  vary  slowing  down  or  acceleration  of  the  flow). 

To  calculate  the  flows  of  this  type  (Figures  la  and  lb),  we  are  going  to  develop  herein  a 
simplified  method  by  using  a  new  variant  of  the  phenomenologic  semi-empiric  turbulent 
boundary-layer  theory  (Novozilov  [1])  that  is  founded  on  the  analogy  with  the  rheologic  power- 
laws  largely  utilized  in  the  study  of  non-Newtonian  non-linear  flows  and  on  the  use  of  the 
Karman  turbulence  model.  It  is  to  note  also  that  the  viscous  sublayer  is  neglected  as  well  as  the 
layer  situated  over  the  zone  for  which  the  universal  logarithmic  law  is  usually  used. 


2.  Planar  (or  axisymmetric)  turbulent  boundary  layer 

The  equations  of  the  planar  turbulent  boundary  layer  can  be  written  in  the  following 
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where  (x,  y)  -  usual  orthogonal  curved  coordinates,  (u,  v)  -velocity  components  at  (x,  y)  in  the 
boundary  layer,  ue(x)-  external  velocity  out  of  the  frontier  of  the  boundary  layer, 

xt  =  |i  —  -pu'v'  =  |i  —  +  t-  global  turbulent  stress  where  the  viscous  part  is  neglected. 

oy  oy 


Fig.  1.  Model  of  divergent  rectangular  channel  (a),  axisymmetrical  diffusor  (b) 
and  meridian  line  shapes  (c). 
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Due  to  the  rheological  power  law  and  the  Prandtl  turbulence  model  chosen  herein,  it 
results  that: 


1  ,  Tn  ou 

—  T  =  vknT  —  , 
P  5y 
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cu 
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dy 

The  system  of  equations  (1)  must  be  resolved  for  the  following  boundary  conditions: 
_  _  5u 

u  =  v  =  0, - >oo,  for  y  =  0 

dy 

u  =  ue(x),  t  =  0,  for  y  =  5(x) 


(2) 

(3) 

(4) 


as  well  as  for  the  next  initial  condition: 

u  =  u0(y)  for  x  =  x0,  (5) 

where  5(x)  is  the  thickness  of  the  boundary  layer. 

In  order  to  calculate  the  turbulent  boundary  layer,  we  used  in  practice  most  often  an 
integral  relation  obtained  from  the  system  of  equations  (1),  called  the  momentum  equation,  in 
the  following  form: 


d52  u'S2/,  it>  1 

5  —  (2  +  H)  =  —  cf 
dx  up  2 


where 


51(x)=f(l-—  )dy,  52(x)=[— (1-  —  )dy,  c 

J  11  J  11  11 


0  Ue 


2T0  rr  _  ^ 

f  “  T’  H  “  s  ' 

pue  °2 


(6) 
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By  analogy  with  the  Novozilov’s  procedure  [1],  the  practical  approximate  method  of  com¬ 
puting  the  turbulent  boundary  layer  proposed  herein  consists  in  the  following  steps  to  be  done: 

-  At  first,  for  a  given  distribution  of  the  external  velocity  ue(x) ,  a  Riccati  like  differential 
equation: 


dz 


v 


,  u: 
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- —  23, 1 63  —  (U-  y  z-  +  2, 097  — z  =  0, 063  ^ 
dx  u  l  ue  v 

is  to  be  integrated  with  a  prescribed  initial  condition  of  52  for  x  =  Xj  : 

z(x)  =  (Ue^2  )4/3  =  zx  for  x  =  Xj. 
v 

-  After  finding  z(x),  determine  the  momentum  thickness  by  using  (10): 

52(x)  =  — z3/4, 


(8) 


(9) 


(10) 


as  well  as  the  parameter  Q(x),  given  by: 


Q(x)  =  —  u’z(x) 


(11) 


-  The  next  step  is  to  determine  the  functions  E(Q),  G(Q)  and  H(Q),  defined  as  follows: 
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E(Q)  =  (2  — n)j[l  +  H(Q)]Q-  io(Q)j  +  2Q,  (12) 

G(Q)  =  cf(-^y-n,  (13) 

V 

H(Q)=|L.  (14) 

°2 

It  is  to  be  noticed  that  the  function  E(Q)  depends  on  n,  not  only  explicitely,  but  also 
through  functions  G(Q)  and  H(Q),  the  forms  of  which  depend  on  n  and  can  be  determined 
numerically.  In  this  sense,  it  is  possible  to  ascertain  that  for  every  particular  value  of  n,  there  is 
“its  proper  formula’’  (12).  In  the  particular  case  (n  =  2/3,  kn  =0,55)  and  for  the  Prandtl 
turbulence  model  (3),  we  carried  out  detailed  computations  of  the  functions  E(Q),  G(Q)  and 
H(Q)  with  approximate  analytical  expressions  of  the  next  forms: 

E(Q)  =  -0.063  +  4.097Q  -  23. 163Q2 

G(Q)  =  0.0938  +  2.1 143Q  +  36. 1035Q2  +  984.7348Q3  >  (15) 

H(Q)  =  1. 6529-2. 17Q 

-  So,  the  formulas: 

H(x)  =  H(Q),  cf(x)  =  z-1/4G(Q)  (16) 

offer  two  particularly  important  characteristics  of  the  turbulent  boundary  layer. 

-  Finally,  by  using  the  well  known  simplest  idea,  it  will  be  possible  to  calculate  an 
improved  approximation  of  the  mean  velocity  in  the  turbulent  boundary  layer  by  ( 1 7): 

-  r  u  ,  io/2)(h-d 
u  _  y  h-1 


ue  L52  H(H  +  1)J 

where  the  functions  52(x)  and  H(x)  are  determined  by  the  formulas  (10)  and  (16). 

Applying  the  analogous  procedure  as  above  in  the  planar  case,  we  can  obtain  also  the 
practical  approximate  method  for  the  axially  symmetric  turbulent  boundary  layer.  The  only 

4  r' 

difference  appears  in  Eq.  8  where  a  new  term - z  should  be  added  at  left  side  (r  is  the  radius 

3  r 

of  the  sections  taken  at  right  angles  to  the  axes  of  revolution). 

3.  Elementary  case 

3.1.  Planar  turbulent  boundary  layer  between  two  divergent  surfaces.  It  is  the  question 
of  two  divergent  surfaces  (i.e.  two  radial  plates  for  k  =  1,  Fig.  la)  in  the  flow  produced  by  a 
linear  source,  so  that  the  external  velocity  can  be  written  in  the  form: 

ue(x)  =  U0(l  +  -Vc,  x  >  0,  k  >  0.  (18) 


If  we  replace  (18)  into  (8),  then  it  results: 

^  =  23, 1 63  (1  +  f  )k-2  z2  +  2, 097  ^  (1  +  f )-'  z  +  0, 063  ^  (1  +  ^)-k , 

dx  LRP  L  L  L  L  L 
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where  Re  =U0L/v.  First,  we  found  the  analytical  solution  of  the  Riccati  differential  equation 
(19),  for  every  other  values  of  k,  except  for: 

0,18139<k<l, 46844,  (20) 


in  the  form: 

T  R  x 

z(x)  = - \y(l  +  -)2-k  M(x)  +  N(x), 

23,163k"  L 


(21) 


M(x)  =  i-  1 


(w  —  1)  —  Cw(w  +  1)(L  +  x)v 


1-U48^  x,_k 

23,163k2  e  L 


2  L  +  x  l  +  Cw(L  +  x)w 
where  the  integration  constant  C,  in  accordance  with  the  initial  condition  z(0)  =  0,  is  the  following: 


C=  W  3-Q97k  +  1  L~w,  with  w  =  (1-6, 194k  +  3, 75432k2)1/2. 
w(w  +  3,097k-l) 


Then,  for  the  values  of  k,  situated  into  the  interval  (20),  we  found  the  solution  in  the  form: 

(22) 

0,06685k-0, 02158 


Q(x)  =  0, 06685  -  0,02158  -  ptg 
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3.2.  Turbulent  boundary  layer  in  an  axisymmetric  diffusor.  In  reality,  it  is  the  question 
of  the  flow  produced  by  a  linear  source  on  a  disc  in  the  presence  of  a  perpendicularly  placed 
funnel  (Fig.  lb).  The  equation  of  the  meridian  line  of  the  funnel  is: 

y=y0(i+|-)k~1,  (23) 

K 

where  y0  =  y(0),  x  >  0,  while  the  distributions  of  ue(x)  and  r(x),  for  different  values  of  the 
exponent  k  (Fig.  lc),  are: 

uc(x)  =  U0(l+-V\  r(x)=R  +  x  ,  (24) 

K 

where  U0  =  ue(0)  and  x  >  0. 

Repeating  the  analogous  procedure  as  in  §3.1,  we  determined  the  both  functions  z(x) 
and  Q(x). 


4.  Calculation  of  the  global  turbulent  boundary  layer  characteristics 

After  determining  the  functions  z(x)  and  Q(x),  we  can  calculate  now  all  global  boundary- 
layer  characteristics  in  accordance  with  the  order  of  operations  announced  at  §2.  So,  in  the  case  of 
the  planar  diffiisors,  we  are  choosing  herein  k  =  1,  i.e.  the  case  of  two  radial  divergent  plates 
(Fig.  la).  Or,  while  this  value  k  =  1  is  situated  over  the  interval  (20),  then  the  solution  is  done  by 
(22),  with  p=0,0258  and  C0  =0,17L.  Next,  the  relation  (11)  offers  z(x),  then  the  momentum  thick¬ 
ness  52(x)  is  obtained  from  (10).  Afterwards,  the  calculation  goes  on  in  accordance  with  the  or¬ 
der  of  operations  above-mentioned.  The  turbulent  boundary  layer  characteristics,  calculated  in  this 
way  in  the  particular  case  Re  =  5.105,  are  presented  in  Figs.  2  and  3.  It  will  be  particularly  re¬ 
marked  that,  the  turbulent  boundary  layer  separation  (G(Q)=0)  occurs  nearly  at  x  =  0,465L.  It  is  to 
be  noticed  that  in  the  analogous  laminar  case  [3],  the  separation  point  appeared  at  x  =  0,161L. 
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Fig.  2.  Variation  of  adimensional 
momentum  thickness. 


Fig.  3.  Variation  of  skin- 
friction  adimensional 
coefficient. 
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Fig.  4.  Comparison  of  the 
detachment  point. 


Concerning  the  axisymmetric  diffusors,  we  are  choosing  herein  k=l,  i.e.,  after  (23)  and 
(24),  the  case  of  an  annular  disc  in  the  flow  produced  by  a  linear  source  (Fig.  lb).  The  turbulent 
boundary  layer  characteristics,  calculated  in  a  similar  way  as  above,  in  the  particular  case 

Re  =  5.105,  are  given  by  Figs.  2  and  3.  The  turbulent  boundary  layer  separation  (G(Q)  =  0) 
appears  at  x  =  0,6R,  while  in  the  analogous  laminar  case  [3],  the  point  of  boundary  layer 
separation  occurred  at  x  =  0,178R. 


5.  Comparison  with  the  other  methods  and  models  of  turbulence 

Figure  4  shows  our  results  [4]  concerning  the  position  of  detachment  points,  together  with 
the  corresponding  results  obtained  by  two  models  of  turbulence:  CFD-k-s-RNG, 
CFD-k-s  by  Rothe  et  al.  [5],  and  UIM  of  Stanford  University  by  Johnston  [6]  as  well  as  the 
experimental  results  of  Ashjaee  et  al.  [7]  It  can  be  remarked  that  the  approximate  method 
predicts  the  detachment  point  with  the  accuracy  of  6  %  approximately  for  the  angles  10°  to  18°. 


6.  Remarks  on  some  (divergent)  external  flows 

The  simplified  practical  method  is  applicable  for  both  internal  and  external  flow 
problems.  So,  an  approximate  procedure  is  proposed  to  evaluate  the  evolution  of  the  velocity 
profile  in  a  neighborhood  of  the  separation  point  of  the  turbulent  boundary  layer  in  the  case  of 
the  Prandtl  turbulence  model  [8].  The  Fig.  5  gives  the  calculated  velocity  profile  for  one  section 
after  the  separation  point  compared  to  the  experimental  results  (presented  by  the  black  points) 
of  Mises  in  the  case  of  the  ident  3800  [2],  while  the  Fig.  6  shows  the  calculated  velocity  profile 
by  the  developed  simplified  method  quite  nearly  the  wall  in  the  same  case. 

7.  Conclusion 

By  beginning  with  the  integral  momentum  equation  and  with  the  analogy  with  the 
rheological  power  law,  which  is  effectivelly  used  to  study  flow  of  non-Newtonian  liquids  with 
nonlinear  viscosities,  a  simple  practical  method  is  obtained  for  an  approximate  computation  of 
both  turbulent  boundary  layers,  planar  and  axisymmetric,  in  the  case  of  the  Prandtl  model  of 
turbulence,  where  for  the  process  ending  it  is  necessary  to  know  only  two  empirical  constants 
(for  example,  n  =  2/3  and  kn  =0,55).  The  developed  method  is  applicable  to  some  internal 
flows,  such  as  a  divergent  rectangular  channel  (Fig.  la)  or  an  axisymmetrical  diffusor  (Fig.  lb), 
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Fig.  5.  Calculated  velocity 
profile  compared  to  the 
experimental  results. 


Fig.  6.  Calculated  velocity 
profile  nearly  to  the  wall. 


Fig.  7.  Turbulent  boundary 
layer  separation  on  the  side 
wall  and  on  the  bottom  of 
the  channel. 


as  well  as  other  planar  or  axisymmetrical,  rectilinear  or  curvilinear,  diffusors  [4].  So,  the  Fig.  7 
summarizes  some  obtained  results  [8]  in  the  case  of  the  divergent  rectangular  channel  for 
Re  =  5.10s.  From  there,  the  turbulent  boundary  layer  separation  appears  earlier  on  the  side  wall 
(dotted  line)  than  on  the  bottom  of  the  channel,  independently  of  the  slowing  down  of  the  flow, 
which  is  in  accordance  with  the  same  tendency,  well-known  before  [3],  in  the  analogous 
laminar  case. 

The  method  gives  relatively  good  results  for  all  global  characteristics  S2 ,  FI  and  cf  of  the 
turbulent  boundary  layer  in  both,  internal  and  external,  problems.  As  expected,  the  velocity 
profiles  are  of  the  lower  quality,  relatively  schematized.  In  return  for  this  drawback,  the  method 
enables  a  satisfactory  evolution  evaluation  of  the  flow  velocity  profile  just  upstream  and 
downstream  the  separation  point  of  the  turbulent  boundary  layer  (Figs.  5  and  6).  Finally,  it  can 
be  noted  that  we  did  also  some  comparaisons  between  the  actual  practical  simplified  method 
and  the  corresponding  experimental  results  [7]  with  satisfactory  concordance. 
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